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Abstract 

We present the first evidence of a retrograde orbit of the transiting exoplanet HAT- 
P-7b. The discovery is based on a measurement of the Rossiter-McLaughlin effect with 
the Subaru HDS during a transit of HAT-P-7b, which occurred on UT 2008 May 30. 
Our best-fit model shows that the spin-orbit alignment angle of this planet is A = 

— 132.6° ti6;3o. The existence of such a retrograde planet have been predicted by 
recent planetary migration models considering planet-planet scattering processes or 
the Kozai migration. Our finding provides an important milestone that supports such 
dynamic migration theories. 
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1. Introduction 

One of the surprising properties of extrasolar planets is their distributions around their 
host stars. Since many Jovian planets have been found in the vicinity (far inside the snow 
line) of their host stars, numbers of theoretical models have been studied to explain inward 
planetary migration. Recently understanding of planetary migration mechanisms has rapidly 
progressed through observations of the Rossiter-McLaughlin effect (hereafter the RM effect: 
Rossiter 1924, McLaughlin 1924) in transiting exoplanetary systems. The RM effect is an 
apparent radial velocity anomaly during planetary transits. By measuring this effect, one 
can learn the sky-projected angle between the stellar spin axis and the planetary orbital axis. 



Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory 
of Japan. 
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denoted by A (see Ohta et al. 2005; Gaudi & Winn 2007 for theoretical discussion). 

So far, spin-orbit alignment angles of about 15 transiting planets have been measured 
(Fabrycky & Winn 2009, and references therein). Among those RM targets, significant spin- 
orbit misalignments have been reported for 3 transiting planets: X0-3b (Hebrard et al. 2008; 
Winn et al. 2009a), HD80606b (Moutou et al. 2009; Pont et al. 2009; Winn et al. 2009b), and 
WASP-14b (Johnson et al. 2009). These misaligned planets are considered to have migrated 
through planet-planet scattering processes (e.g., Rasio & Ford 1996; Marzari & Weidenschilling 
2002; Nagasawa et al. 2008; Chatterjee et al. 2008) or Kozai cycles with tidal evolution (e.g., 
Wu & Murray 2003; Takeda & Rasio 2005; Fabrycky & Tremaine 2007; Wu et al. 2007), rather 
than the standard Type II migration (e.g., Lin & Papaloizou 1985; Lin et al. 1996; Ida & Lin 
2004). 

The existence of such misaligned planets has demonstrated validity of the planetary 
migration models considering planet-planet scattering or the Kozai migration. On the other 
hand, such planetary migration models also predict significant populations of "retrograde" 
planets. Thus discoveries of retrograde planets would be an important milestone for confirming 
the predictions of recent planetary migration models, and intrinsically interesting. 

In this letter, we report the first evidence of such a retrograde planet in the transiting 
exoplanetary system HAT-P-7. Section 2 summarizes the target and our Subaru observations, 
and section 3 describes the analysis procedures for the RM effect. Section 4 presents results and 
discussion for the derived system parameters. Finally, section 5 summarizes the main findings 
of this letter. 

2. Target and Subaru Observations 

HAT-P-7 is an F6 star at a distance of 320 pc hosting a very hot Jupiter (Pal et al. 2008; 
hereafter P08). Among transiting-planet host stars, F type stars are interesting RM targets 
because these stars often have a large stellar rotational velocity, which facilitates measurements 
of the RM effect. However, the rotational velocity of HAT-P-7 is V^sin = 3.8 km s"^ (P08), 
which is unusually slower than expected for an F6 type star. Nevertheless, this system is 
favorable for the RM observations, since the star is relatively bright {V = 10.5) and the expected 
amplitude of the RM effect {V sin Is{Rp/ Rs)"^ ~ 20 m s~^) is sufficiently detactable with the 
Subaru telescope. 

We observed a full transit of HAT-P-7b with the High Dispersion Spectrograph (HDS: 
Noguchi et al. 2002) on the Subaru 8.2m telescope on UT 2008 May 30. We employed the 
standard I2a set-up of the HDS, covering the wavelength range 4940 A < A < 6180 A and used 
the Iodine gas absorption cell for radial velocity measurements. The slit width of O'.'G yielded 
a spectral resolution of ~60000. The seeing on that night was around 0^'6. The exposure time 
for radial velocity measurements was 6-8 minutes, yielding a typical signal-to-noise ratio (SNR) 
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of approximately 120 per pixel. We processed the observed frames with standard IRAF^ proce- 
dures and extracted one-dimensional spectra. We computed relative radial velocities following 
the algorithm of Butler et al. (1996) and Sato et al. (2002), as described in Narita et al. (2007). 
We estimated the internal error of each radial velocity as the scatter in the radial-velocity so- 



radial velocities and uncertainties are summarized in table 1. 
3. Analyses 

We model the RM effect of HAT-P-7 following the procedure of Winn et al. (2005), as 
described in Narita et al. (2009) and Hirano et al. in prep. We start with a synthetic template 
spectrum, which matches for the stellar property of HAT-P-7 described in P08, using a synthetic 
model by Coelho et al. (2005). To model the disk- integrated spectrum of HAT-P-7, we apply a 
rotational broadening kernel of ^sin = 3.8 km s~^ and assume limb- darkening parameters for 
the spectroscopic band as Ui = 0.45 and U2 = 0.31, based on a model by Claret (2004). We then 
subtract a scaled copy of the original unbroadened spectrum with a velocity shift to simulate 
spectra during a transit. We create numbers of such simulated spectra using different values 
of the scaling factor / and the velocity shift Vp, and compute the apparent radial velocity of 
each spectrum. We thereby determine an empirical formula that describes the radial velocity 
anomaly Av in HAT-P-7 due to the RM effect, and find 



For radial velocity fitting, including the Keplerian motion and the RM effect, we adopt 
stellar and planetary parameters based on P08 as follows; the stellar mass Mg = 1.47 [Mq], the 
stellar radius Rg = 1.84 [i?©], the radius ratio Rp/Rs = 0.0763, the orbital inclination i = 85.7°, 
and the semi-major axis in units of the stellar radius a/Rs = 4.35. We assess possible systematic 
errors due to uncertainties in the fixed parameters in section 4. We also include a stellar jitter 
of 3.8 m s~^ for the P08 Keck data as systematic errors of radial velocities by quadrature sum. 
It enforces the ratio of contribution and the degree of freedom for the Keck data to be unity. 
We do not include additional radial velocity errors for the Subaru data, because we find the 
ratio for the Subaru dataset is already smaller than unity (as described in section 4). 

In addition, we adopt the transit ephemeris Tc = 2453790.2593 [HJD] and the orbital 
period P = 2.2047299 days based on P08. Note that this ephemeris has an uncertainty of 3 
minutes for the observed transit; however the uncertainty is well within our time- resolution 
(exposure time of 6-8 minutes and readout time of 1 minute) and is negligible for our purpose. 
The adopted parameters above are summarized in table 2. 

^ The Image Reduction and Analysis Facility (IRAF) is distributed by the U.S. National Optical Astronomy 
Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 



lutions among ~4 A segments of the spectrum. The typical internal error was ~5 m s 



\ The 
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Our model has 3 free parameters describing the HAT-P-7 system: the radial velocity 
semi-amplitude K, the sky-projected stellar rotational velocity V smls-, and the sky-projected 
angle between the stellar spin axis and the planetary orbital axis A. We fix the eccentricity e 
to zero, and the argument of periastron uj is not considered. Finally we add two offset velocity 
parameters for respective radial velocity datasets {vi. our Subaru dataset, V2'- the Keck dataset 
in P08). 

We then calculate the statistic (hereafter "main-case") 

n 2 



ai J 

where Wj^obs and cTj are observed radial velocities and uncertainties, and f i,caic are radial velocity 
values calculated based on a Keplerian motion and on the empirical RM formula given above. 

We determine optimal orbital parameters by minimizing the statistic using the 
AMOEBA algorithm (Press et al. 1992). We estimate la uncertainty of each free parame- 
ter based on the criterion Ax^ = 1.0 when a parameter is stepped away from the best-fit value 
and the other parameters are re-optimized. 

We also fit the radial velocities using another statistic function for reference (hereafter 
"test-case"), 



n 2 



l^sinL-3.8 
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(3) 



0.5 

The last term is a priori constraint for V^sin Jg to match the independent spectroscopic analysis 
by P08. 

4. Results and Discussion 

Figure 1 shows observed radial velocities and the best-fit model curve for the main- 
case. Figure 2 illustrates the RM effect of HAT-P-7b with the best-fit model and also shows 
a comparison with the case of A = 0° and Vsinlg = 3.8 km s~^ model. The upper panel 
of figure 3 plots a contour in {X,V sinlg) space. As a result, we find the key parameter 
A = —132.6° (+10.5°, —16.3°), implying a retrograde orbit of HAT-P-7b. The stellar rotational 
velocity is Vsinlg = 2.3 (+0.6,-0.5) km s~^, which is marginally consistent with the P08 
spectroscopic result (V^sin/g = 3.8 + 0.5 km s~^). Residuals from the best-fit model indicate 
rms of 4.14 m s~^ for the Subaru dataset and 4.09 m s~^ for the P08 Keck dataset. The rms 
of the Subaru residuals is well within our internal radial velocity errors, and that of the Keck 
residuals is in good agreement with the assumed jitter level of 3.8 m s~^. One may wonder 
that a smaller Vsinlg allowed in the main-case would weaken the detection-significance of the 
RM effect. However, since Vsinls = km s~^ is excluded by Ax^ = 36.49, there is very little 
chance that a true V sin Is is actually nearly zero and a spin-orbit alignment angle A is very 
small. The lower panel of figure 3 plots a similar contour but for the test-case. In this case, 
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we find A = —122.5° (+6.4°, —7.7°) and ^ sin = 3.1 ± 0.4 km s~^. Thus our two results (main 
and test cases) are well consistent with each other. In addition, we test the fitting with the 
eccentricity e and the argument of periastron uj as free parameters. As a result, we do not find 
any significant (nonzero) eccentricity for this planet. The best-fit parameters and uncertainties 
are summarized in table 3. 

In the above analyses, we fixed several parameters as summarized in table 2, which were 
based on P08 and Claret (2004). In order to estimate the level of possible systematic errors, we 
retry the fitting for following four cases; (1) a/ Rs = 4.63, i = 89.2° (corresponding to la lower 
limit of the impact parameter in P08); (2) a/Rg = 3.97, i = 82.6° (corresponding to la upper 
limit of the impact parameter in P08); (3) ui = 0.65 (a greater limb-darkening case); and (4) 
Ml = 0.25 (a smaller limb- darkening case). 

Consequently, we find that respective results for A and sin are; (1) A = 
-151.8° (+11.1°, -12.8°) and V sin J, = 2.1 ± 0.4 km s"^; (2) A = -99.5° (+3.0°, -6.2°) and 
sin/, = 7.6 + 2.9 km s'^ (3) A = -135.9° (+11.2°, -16.7°) and 1/ sin/, = 2.2 + 0.5 km s^^ and 
(4) A = -129.8° (+9.9°, -15.6°) and Vsinls = 2.4 + 0.6 km s'^ Thus A = -164.6° - -96.5° 
and V sin Is = 1.7 - 10.5 km s~^ can be still probable if the uncertainties for fixed parameters 
(especially for the impact parameter) are taken into account. These systematic errors would 
be significantly reduced when the Kepler photometric data for HAT-P-7 are available (Borucki 
et al. 2009). 

The derived value of A seems to indicate a retrograde orbit by itself. However, since the 
true spin-orbit angle |\E'| also depends on the inclination of the stellar spin axis /,, |\&| is not 
necessarily greater than 90° (corresponding to a retrograde orbit) even if A = —132.6°. Thus 
one might wonder whether the planet HAT-P-7b is statistically in a retrograde orbit. We can 
roughly estimate the probability using the relation of spherical geometry, 

cos I \I/ 1 = COS Is COS i + sin /, sin i cos A. (4) 

Note that /, ranges from 0° to 180°. We compute the true spin-orbit angle of HAT-P-7b 
by substituting the observed values of i and A into the relation. We adopt i = 85.7° (P08) 
and test three representative cases for A (—164.6°, —132.6°, and —96.5°). Assuming an uniform 
distribution for cos/, within the range of value, the probabilities of a retrograde orbit (|\&| > 90°) 
are 99.85% (A = -164.6°), 99.70% (A = -132.6°), and 91.65% (A = -96.5°), respectively. We 
note that is always larger than 85.7° (the adopted value of i, in the case of /, = 0°). Those 
estimates favor a retrograde orbit of HAT-P-7b. 

On the other hand, it is important to point out that the stellar rotational velocity 
V sin Is = 3.8 km s^^ determined by the spectroscopic analysis (P08) is exceptionally slow 
for an F6 star. For example, HAT-P-2 and TrES-4, which are other known planetary host 
stars with similar spectral type, have larger stellar rotational velocities: V sin Is = 19.8 km s~^ 
(HAT-P-2, Bakos et al. 2007 and also confirmed from the RM effect by Winn et al. 2007) and 
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V sin Is = 8.5 km (TrES-4, Sozzetti et al. 2009 and also confirmed from tlie RM effect by 
Narita et al. in prep.). The small V^sinJ^ may suggest a smaller inclination angle of the stellar 
spin axis. In that case, it is highly possible that the planet is in a nearly polar retrograde orbit. 

We note that too small Ig can be constrained by the facts that a faster stellar rotation 
of HAT-P-7 than 500 km s~^ would be physically unlikely due to stellar break-up, and that 
a faster rotation than 100 km s~^ would be empirically unlikely for a F6 star (Gray 2005). 
Translating the constraints into cos/^, we find a probability of such unrealistic cases is only 
~0.03%, which has very little impact on the probability estimations for a retrograde orbit. In 
any case, it would be important to directly constrain Ig by other observational methods (e.g., 
Henry & Winn 2008 or asteroseismology with the Kepler) in order to estimate a true spin-orbit 
angle of HAT-P-7b. 

We previously experienced a false positive of a spin-orbit misalignment in HD17156b due 
to lower precision radial velocity data (Narita et al. 2008; Narita et al. 2009). The problem for 
the HD17156 case (Narita et al. 2008) was that radial velocity uncertainties were comparable 
with predicted RM amplitude and also due to a poor number of radial velocity samples. Based 
on the lesson, we estimate the significance of our RM detection using the equation (26) in Gaudi 
& Winn (2007). As a result, the SNR of our RM detection is over 10, and thus we conclude to 
have obtained radial velocities of a sufficient number and precision to model the RM effect of 
HAT-P-7b. 

Nevertheless, we finally note that we should care about possible systematic errors in A. 
Since the RM amplitude of HAT-P-7b is only ~ 15 m s~^, any systematic shift as much as 
several m s~^ due to stellar jitter or other reasons at ingress or egress phase would cause a 
large systematic difference in A. Thus further radial velocity measurements for this interesting 
system are desired to confirm a retrograde orbit of HAT-P-7b more decisively. 

5. Summary 

We observed a full transit of HAT-P-7b with the Subaru 8.2m telescope on UT 2008 
May 30, and measured the RM effect of this planet. Based on the RM modeling, we discovered 
the first evidence of a retrograde orbit of HAT-P-7b. This is the first discovery of a retrograde 
extrasolar planet. The existence of such planets have been indeed predicted in some recent 
planetary migration models considering planet-planet scattering and/or the Kozai migration 
(e.g., Fabrycky & Tremaine 2007; Wu et al. 2007; Nagasawa et al. 2008; Chatterjee et al. 2008). 
In addition, it is interesting to point out that HAT-P-7b is the first spin-orbit misaligned planet 
having no significant eccentricity. This discovery may suggest that other hot Jupiters in circular 
orbits also have significant spin-orbit misalignments or even retrograde orbits. Thus further 
RM observations for transiting planets, including not only eccentric or binary system planets 
but also close-in circular planets, would be encouraged in order to understand populations of 
aligned / misaligned / retrograde planets. 
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Note added after submission — Winn et al. (2009c) report an independent evidence for 
a retrograde orbit of HAT-P-7b, based on independent Subaru HDS observations conducted in 
2009 June/July. 

We acknowledge the invaluable support for our Subaru observations by Akito Tajitsu, 
a support scientist for the Subaru HDS. We are grateful to Yasushi Suto, Ed Turner, Wako 
Aoki, and Torn Yamada for helpful discussions on this study; Josh Winn and his colleagues 
for sharing their information in advance of publication. This paper is based on data collected 
at Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. 
The data analysis was in part carried out on common use data analysis computer system at 
the Astronomy Data Center, ADC, of the National Astronomical Observatory of Japan. N.N. 
is supported by a Japan Society for Promotion of Science (JSPS) Fellowship for Research (PD: 
20-8141). We wish to acknowledge the very significant cultural role and reverence that the 
summit of Mauna Kea has always had within the indigenous Hawaiian community. 
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Fig. 1. Upper panels: Radial velocities and the best-fit curve of HAT-P-7 as a function of orbital phase. 
The upper figure show the entire orbit and the lower figure do the zoom of transit phase. Bottom panels: 
Residuals from the best-fit curve. 
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Fig. 2. Upper figure: The RM effect of HAT-P-7b. The upper panel shows difference radial velocities 
(namely, the Keplerian motion is subtracted from the original radial velocities). The solid line indicates 
the best-fit RM model. The lower panel plots residuals from the best-fit model. Lower figure: The case 
for A = 0° and Vsmis = 3.8 km s^^ model, for reference. 
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Table 1. Radial velocities obtained with the Subaru/HDS. 
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Table 1. (Continued.) 
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-75.20 


5.14 



Table 2. Adopted stellar and planetary parameters. 



Parameter 


Value 


Source 


M, [Mq] 


1.47 


P08 


Rs [Rq] 


1.84 


P08 


Rp/Rs 


0.0763 


P08 


^n 


85.7 


P08 


a/Rs 


4.35 


P08 


Ui 


0.45 


Claret (2004) 


U2 


0.31 


Claret (2004) 


Tc [HJD] 


2453790.2593 


P08 


P [days] 


2.2047299 


P08 



Table 3. Best-fit values and uncertainties of the free parameters. 





main-case 


test-case 


Parameter 


Value 


Uncertainty 


Value 


Uncertainty 


K [m s-i] 


212.6 


±1.9 


213.3 


±1.9 


Vsinlg [km s~^] 


2.3 


-1-0.6 
-0.5 


3.1 


±0.4 


A [T 


-132.6 


-1-10.5 
-16.3 


-122.5 


+6.4 
-7.7 


Vi [m s^-'-] 


-14.7 


±1.6 


-16.6 


±1.3 


rms (Subaru) [m s~^] 


4.14 




4.32 




V2 [m s~^] 


-37.4 


±1.6 


-37.5 


±1.6 


rms (Keck) [m s^-'^] 


4.09 




4.09 





~: Systematic errors are not included in the uncertainties (see 
text). 
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